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ABSTRACT 

An  echo  sounder  can  be  used  to  indicate  the  type  of  sediment  on  the  sea  floor  through 
analysis  of  flie  returned  echo's  shape  and  energy  characteristics.  This  document 
describes  the  ABCS  program  written  to  analyse  echosounder  ping  returns  for  simple 
shape  and  energy  statistics,  with  the  intention  to  investigate  which  statistics  are 
important  in  characterising  the  seabed.  ABCS  was  wntten  in  Visual  Basic  and  can  read 
in  the  data  produced  by  the  Edio  Listener  device  manufactured  by  Sonar  Data 
Tasmania.  Also  discussed  in  this  paper  is  the  need  to  transform  all  the  data  obtained  to 
a  reference  seafloor  depth  and  the  method  used  by  ABCS  to  do  this. 
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The  ABCS  Program  for  the  Analysis  of  Echo 
Sounder  Returns  for  Acoustic  Bottom 
Classification 


Executive  Summary 

The  type  of  sediment  found  in  a  shaUow  water  environment  (less  than  200  m  depth) 
has  a  major  effect  on  the  propagation  of  sound  in  the  water  through  repeated  bottom 
interaction.  If  the  propagation  loss  or  acoustic  backscatter  in  an  area  are  to  be  predicted 
the  sediment  types  in  this  area  must  be  known.  Since  the  sediment  type  can  change 
dramaticahy  over  short  distances  in  shallow  water,  the  costs  involved  in  taking  enough 
bottom  samples  and  analysing  these,  even  for  a  small  area,  could  be  massive.  Therefore 
investigating  a  cheaper  method  of  obtaining  the  sediment  type  would  be  worth  while. 

Since  the  returned  signal  from  a  boat's  echo  sounder  is  dependent  on  the  type  of 
sediment,  analysing  the  shape  and  energy  of  the  returned  signal  could  give  the 
sediment  type.  This  would  be  a  far  cheaper  method  to  obtain  the  sediment  type  data. 
Due  to  the  high  frequency  ranges  used  by  most  echo  sounders  this  method  will  usually 
give  the  sediment  type  for  only  the  thin  surface  layer  of  the  sea  floor. 

Analysing  echo  sotmder  returns  has  been  done  before  with  varied  results,  but 
companies  now  doing  this  work  regard  their  processing  as  proprietary,  and  methods 
vary  from  one  manufacturer  to  another. 

To  examine  how  an  echo  sounder  return  changes  with  different  sediment  t5q)es  an 
Echo  Listener™  device  was  attached  to  the  Furuno™  echo  sounder  on  the  DSTO 
Sydney  work  boat  and  die  program  ABCS  was  written  to  analyse  the  returned  wave 
forms. 

ABCS  reads  the  files  produced  by  die  Echo  Listener,  finds  the  seafloor  echo  sounder 
returns  in  the  signal,  and  calculates  a  number  of  statistics  about  them. 
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1.  Introduction 

The  type  of  sediment  found  in  a  shallow  water  environment  (less  than  200  m 
depth)  has  a  major  effect  on  the  propagation  of  soxmd  in  the  water  through  repeated 
bottom  interaction  (D.  Cato  et  al  1995).  If  the  propagation  loss  or  acoustic  backscatter  m 
an  area  are  to  be  predicted  the  sediment  types  in  this  area  must  be  known.  Smce  the 
sediment  type  can  change  dramatically  over  short  distances  in  shaUow  water,  the  costs 
involved  in  taking  enough  bottom  samples  and  analysing  them  even  for  a  small  area 
could  be  massive.  Therefore  investigating  another  method  of  obtaining  the  sediment 
type  which  is  cheaper  would  be  worth  while. 

Since  tire  returned  signal  from  a  boat's  echo  sounder  is  dependent  on  the  type  of 
sediment,  analysing  the  shape  and  energy  of  the  returned  signal  could  give  the 
sediment  type.  This  would  be  a  far  cheaper  method  to  obtain  the  sediment  t5^e  data. 
Due  to  the  high  frequency  ranges  used  by  most  echo  sounders  this  method  will  usually 
give  the  sediment  type  for  a  thin  surface  layer  of  the  seafloor,  similar  to  the  results 
found  from  a  sediment  grab  sample.  If  the  sediment  types  were  required  deeper  into 
the  seafloor  then  a  core  sample,  or  an  inversion  technique  would  have  to  be  used. 

Analysing  echo  soimder  returns  has  been  done  before,  with  varied  results.  The 
main  commercial  products  available  now  are  QTC-View™  and  RoxAnn™.  RoxAnn 
uses  die  energy  from  the  first  and  second  returns^  to  give  the  sediment  type  (Bums  et 
al  1989),  while  QTC-View  uses  a  large  number  of  echo  shape  parameters  from  the  first 
return  only  to  give  the  sediment  type  (Prager  et  al  1995).  To  obtain  the  sediment  type 
both  require  calibration,  which  involves  transiting  over  known  sediment  types  at  a 
similar  depth  to  the  area  that  needs  to  be  surveyed. 

The  problem  with  these  two  products  is  that  the  companies  now  doing  this  work 
regard  their  processing  as  proprietary.  They  give  out  little  information  on  what 
calculations  tiie  units  are  doing  to  get  a  sechment  type  (this  is  done  to  protect  the 
companies'  patents).  Also  there  is  a  high  original  purchase  cost,  and  the  reliance  on  the 
companies  to  keep  upgrading  the  product  along  with  the  cost  of  upgrades.  The  current 
cost  of  the  RoxAnn  system  is  around  $20,000  per  unit,  with  QTC-View  costing  between 
$20,000  to  $100,000  per  unit  depending  on  what  software  is  purchased.  Since  little 
information  is  given  out  on  what  calcxdations  are  being  done  to  get  the  sediment  types, 
the  theory  behind  the  calculations  can  not  be  investigated  and  only  the  manufacturers 
can  do  upgrading  of  the  detection  methods. 

To  try  to  overcome  these  problems  and  improve  DSTO's  rmderstanding  of  the 
mathematics  used  to  obtain  the  sediment  types,  an  Edio  Listener  device  manufactured 
by  Sonar  Data  Tasmania  Pty  Ltd  was  attached  to  the  Furuno™  echo  sounder  on  the 

1  The  first  return  has  a  single  bottom  bounce  (from  a  source  /  seafloor  /  receiver  path),  while 
the  second  return  has  two  bottom  bounces  (from  a  source  /  seafloor  /  sea  surface  /  seafloor  / 
receiver  path). 
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DSTO  Sydney  work  boat  to  record  die  echo  waveforms.  A  program  (called  ABCS.exe) 
was  tiien  written  to  analyse  the  returned  waveforms.  This  report  describes  the  running 
of  the  ABCS  program,  what  shape  statistics  it  produces,  and  what  needs  to  be 
considered  when  running  the  program.  The  current  version  of  the  ABCS  program  is 
not  trying  to  obtain  the  sediment  type.  It  was  written  to  see  how  selected  shape  and 
energy  statistics  from  the  returned  echo  sounder  signal  change  with  sedirnent  type, 
thus  determining  which  shape  statistics  are  important  for  sediment  classification. 


2.  The  Data 


The  data  used  to  test  the  ABCS  program  was  produced  using  an  Echo  Listener™ 
Device.  This  data  was  suppUed  in  CSV  (comma  separated  variable)  format,  and  was 
obtained  at  a  sampling  frequency  of  13.6  kHz.  The  echo  sounder  depth  setting  was  set 
to  110  m  giving  2000  samples  per  ping.  If  the  data  format  is  changed  from  CSV  a  slight 
modification  will  have  to  be  made  to  the  ABCS  program. 

The  Echo  Listener  device  is  an  A/D  converter  that  has  been  connected  in  parallel  to 
the  Furuno™  echo  sounder  on  the  DSTO  Sydney  work  boat.  The  Echo  Listener  device 
can  continuously  log  the  returned  echo  sounder  signal  onto  a  PC  or  laptop  without 
affecting  the  operation  of  the  echo  sounder.  The  Echo  Listener  was  manufactured  and 
supplied  by  Sonar  Data  Tasmania  Pty  Ltd^ 


3.  The  'ABCS'  Program 


The  ABCS  program  was  written  in  Visual  Basic  and  consists  of  five  display  forms. 
These  are  the  main  display  form  (which  is  used  to  enter  the  data),  three  plot  displays, 
and  a  statistics  display  form.  Only  the  main  display  form  needs  to  be  open,  the  others 
can  be  opened  and  closed  as  required. 

This  program  was  not  designed  to  process  the  pings  in  real  time,  but  reads  the 
output  files  from  the  Echo  Listener™  device  and  calculates  the  statistics  from  them  in 
post  processing.  The  current  version  will  be  used  to  examine  how  selected  shape  and 
energy  statistics  change  with  sediment  type;  it  is  not  trying  to  classify  the  sediment 
types. 


2  For  more  information  on  Sonar  Data  and  tiie  Echo  Listener  see  their  web  site  at 
http:/  /  www.verdant.com.au/SonarData/ 
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3.1  Depth  Nonnalisation 

The  shape  and  power  of  the  returned  signal  from  the  seafloor  can  change 
significantly  with  seafloor  depth  as  a  result  of  spreading  losses  ^d  absorption  (see 
appendix  1),  even  if  the  seafloor  sediment  remains  the  same.  Seafloor  dep&  ch^ges 
also  cause  dilation  or  compression  (see  appendix  1)  of  the  returned  pmg.  ese 
changes  due  to  seafloor  depth  variation  will  effect  the  shape  statistics  bemg  calculated 
by  flie  ABCS  program,  thus  obscuring  the  changes  due  to  variations  m  sediment  type. 
To  overcome  this,  the  returned  signal  is  transformed  to  a  reference  seafloor  depth  and 
the  statistics  calculated  at  the  reference  seafloor  depth. 

To  transform  a  returned  signal  to  the  reference  seafloor  depth,  time  ^d  power 
adjustments  need  to  be  made.  The  time  adjustment  is  first  made  to  adji^t  the  duration 
of  the  returned  ping,  which  changes  due  to  the  differei^  path  lengths  ^ol^^d  m 
backscattering  the  ping  from  different  seafloor  depths.  The  power  adjustment  then 
removes  tire  effect  of  spherical  spreading. 


3.1.1  Time  Adjustment 

The  time  adjustment  enables  returns  from  the  actual  seafloor  (dep&  d)  and  the 
reference  seafloor  (depth  do)  to  maintain  the  same  time/angle  relationships  (Caugfiey 
et  al  1994)  Sampling  at  the  same  angles  for  different  seafloor  depths  removes  the  need 
to  allow  for  beam  patterns,  and  for  the  bottom  backscatter  function  changing  with 
incidence  angle. 


The  time  adjustment  is  ^  ^  (1)  (Caughey  et  al.  1994) 

where  d  =  The  actual  seafloor  depth. 

do  =  The  reference  seafloor  depth. 

Therefore  t'=—=^  (2) 

y  a 

where  f '  =  The  time  adjusted  to  the  reference  depth. 

t  =  The  actual  time  of  the  echo  sounder  signal  firom  the  actual  seafloor  depth. 

These  tirnps  start  at  zero  when  the  echo  sotmder  sends  the  ping,  with  the  time 
adjustment  having  the  effect  of  expanding  or  compressing  the  returned  signal  along 
the  time  axis. 
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3.1.2  Power  Adjustment 

Since  the  backscattered  energy  received  at  the  hydrophone  from  the  seafloor  is 
bottom  reverberation,  this  energy  received  can  be  found  using: 

RL  =  SL-  401ogd  +  5,  +  lOlog^  (D.  Cato  et  al  1995) 

where  RL  =  The  energy  flux  from  the  actual  seafloor  in  dB  re  J/  m^. 

SL  =  The  source  level  energy  of  the  echo  soimder  in  dB  re  J/  sr. 
d  =  The  actual  seafloor  depth  in  m. 

Sb  =  The  boundary  scattering  strength  in  dB. 

A  =  The  area  of  the  insonified  section  of  seafloor  at  the  actual  seafloor  depth 
inm^. 

A  =  7id^isaa  (see  Figure  1) 

The  above  equation  can  also  be  used  at  the  reference  depth. 

RZ,’=5Z-401ogd’+5'4  +101og^' 

where  RL'  =  The  energy  flux  from  the  reference  seafloor  in  dB  re  J/ m2. 
d'  —  The  reference  seafloor  depth. 

A '  =  The  area  of  the  insonified  section  of  seafloor  at  the  reference  depth. 

A'=  nd'^  tana  (see  Figure  1) 

Sb  and  SL  will  be  the  same  at  both  seafloor  depths. 

By  combining  these  two  equations  the  change  in  energy  flux  with  seafloor  depth 
can  be  foimd. 

RL'-RL  =  40  log  d  -  40  log  d'-l  0  log  ^  + 10  log  ^'  =  20  log 

This  energy  flux  is  converted  into  a  voltage  by  the  hydrophone.  Since  the 
hydrophone  size  is  not  changing,  the  area  component  of  the  energy  flux  is  constant  at 
both  depths.  Therefore  the  change  in  energy  can  be  written  as: 


where  E  =  The  received  energy  at  the  hydrophone  from  the  actual  seafloor  in  Joules. 
E'  =  The  adjusted  received  energy  for  the  reference  seafloor  in  Joules. 
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Figure  1:  Geometry  of  the  Pulse  Insonifying  the  Seafloor 

The  Echo  Listener  produces  output  files  containing  the  power  levels  returned  from 
the  seafloor,  and  since  equation  (2)  gives  a  change  in  the  time  scale,  this  must  be 
considered  when  calculating  the  adjusted  depth. 


dE'  „  dE 
df  dt 

where  1?'=  The  adjusted  received  power. 

R  =  The  actual  received  power. 

Combining  these  two  equations  and  substituting  in  equation  (3)  gives: 
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Since  —  =  r—  (from  equation  (2))  the  above  equation  can  be  written  as: 

dt  dt 


R  ' 
But 


dt 


\^oJ 


d 

vE  — 
'  dt 


r  ^  ^  V 
+rl 


d„ 


d. 


o) 


dt 


dE 

dt 


dE 

dt 


dt 


\^oJ 


=  0  ,  so  the  power  adjustment  required  is: 


R 

R=’' 


f 


(4) 


For  the  first  return  the  transmitted  and  received  angles  are  the  same  (ie.  a),  but  for 
the  second  return  the  received  angle  is  not  the  same  as  the  transmit  angle  (Heald  and 
Pace  1996).  This  does  not  affect  the  power  adjustment  equation,  since  the  time 
adjustment  enables  returns  from  the  actual  seafloor  depth  and  the  reference  seafloor 
depth  to  maintain  the  same  time/  angle  relationship  (Caughey  et  al  1994).  Therefore  the 
time  and  power  corrections  can  be  applied  to  the  whole  returned  signal- 


3.1.3  Effect  of  Transforming  to  a  Reference  Seafloor  Depth 

Simple  adjustments  like  eq(2)  and  eq(4)  do  not  give  an  exact  transformation  to  the 
reference  seafloor  depth.  Even  in  a  simple  case  of  a  flat  seafloor  and  constant 
backscatter  with  angle,  the  above  equations  don't  consider  the  geometric  changes 
exactly  (see  appendix  2). 

If  we  were  to  try  to  precisely  transform  the  echo  sounder  returns  to  a  reference 
seafloor  depth  the  ping  duration,  ping  intensity,  and  Echo  Listener™  sampling  rate 
would  have  to  be  changed  depending  on  the  seafloor  depth.  Since  this  is  difficult, 
simple  transformations  are  done  to  get  the  returned  sigiial  to  a  reference  seafloor 
depth. 

To  test  how  well  the  transformations  in  equations  (2)  and  (4)  work,  a  simple  example 
was  considered.  Using  a  number  of  piece-wise  constant  power  inputs  to  simulate  a  sine 
wave  ping  envelope  (see  Figure  2),  the  program  "Levels"  was  written  to  see  what  the 
returned  signak  would  look  like  (see  appendix  3).  The  program  was  run  for  a  2  msec 
ping  duration  (ttus  duration  was  used  since  it  is  the  ping  duration  of  the  Furuno™  echo 
sounder)  and  the  effect  of  transforming  a  100  m  seafloor  depth  signal  to  a  reference 
seafloor  depth  of  50  m  considered  (see  Figure  3).  This  showed  the  transformed  signal 
had  the  same  arrival  time  and  an  almost  identical  power  variation  with  time  as  the 
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50  m  Signal.  Tlierefore  for  the  Furuno  echo  sounder  this  transformation  should  give  a 
good  result  as  long  as  the  environment  is  close  to  the  assumptions  made  by  Levels. 
These  assumptions  made  by  Levels  were  an  omnidirectional  source  and  receiver,  a  flat 
seafloor,  and  a  backscatter  strength  independent  of  grazing  angle. 


Figure  2:  A  simulated  sine  wave  -ping  made  up  of  a  number  of  piece-wise  constant  segments. 
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Time  (sec) 

Figure  3  :  Returned  Signal  for  a  2  msec  Sine  Wave  Source.  The  100  m  depth  transformed  to  50 
m  and  the  50  m  depth  No  corrections  lines  overlie  each  other.  This  plot  loas 
calculated  using  Levels. 
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If  a  longer  ping  duration  was  used  this  match  would  not  be  as  precise,  but  even  a 
long  echo  sounder  ping  of  200  msec  with  a  large  change  in  seafloor  depth  still  has  a 
good  match  (see  Figure  4). 
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Figure  4:  Returned  Signal  for  a  200  msec  Sine  Wave  Source.  This  plot  was  calculated  using 
Levels. 


3.2  Running  the  ABCS  Program 

The  ABCS  program  needs  the  actual  seafloor  depth,  reference  seafloor  depth,  water 
sormd  speed,  echo  sotmder  maximum  depth  setting,  time  duration  of  the  ping  used, 
and  the  file  containing  the  received  Echo  Listener™  returns.  All  the  data  are  entered  on 
the  Main  form  (Input  Values  frame).  All  the  operations  are  done  on  the  Main  form  (see 
Figure  5),  the  other  forms  are  for  display  only. 
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Figure  5:  The  Main  Form  for  program  ABCS.  The  Actual  Return  Times  and  Reference  Return 
Times  are  calculated  by  ABCS,  these  frames  can  not  be  changed  by  the  operator. 


The  first  step  is  to  enter  fiie  information  in  the  Input  Values  frame,  and  then  open 
the  CSV  file  produced  by  the  Echo  Listener.  The  Sounder  depth  setting  must  be  entered 
before  the  file  is  opened,  since  after  the  file  is  opened  this  value  can  no  longer  be 
changed.  The  other  values  in  the  Input  Values  frame  can  be  altered  during  the  analysis 
as  required.  To  open  the  echo  soimder  file  press  the  Open  button,  and  locate  the 
relevant  file.  The  program  displays  the  file  name  in  the  header  of  the  Main  form.  After 
the  file  name  is  entered  the  program  first  looks  for  the  file  header.  This  is  a  listing  of  all 
the  sampled  depths  for  each  ping.  By  finding  how  many  depth  samples  exist  in  the  file 
header,  the  Echo  Listener  sampling  rate  can  be  found. 


2  X  SounderMaxDepthSetting  .g. 

SamplingRate  x  NumberOfHeaderSamples 

The  numerator  is  2  x  Sounder  Max  Depth  Setting  since  this  is  the  return  path  length. 

The  file  header  is  also  used  to  determine  the  number  of  samples  per  ping,  since 
each  ping  in  a  particular  file  has  the  same  number  of  samples.  Each  ping's  set  of 
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returns  is  analysed  one  at  a  time,  there  is  presently  no  averaging  of  a  number  of  returns 
in  this  version  of  the  program.  When  recording  the  data  the  echo  sounder  smaxirnum 
denth  should  be  set  deep  enough  so  that  the  first  and  second  returns  can  be  recorded 
before  the  next  ping  is  sL  (ie.  the  Sounder  Depth  Setting  should  be  greater  than  twice 

the  seafloor  depth). 

The  program  wUl  then  read  in  the  returns  for  the  first  ping.  The  next 
the  start  and  end  of  the  first  and  second  returns  (The  first  return  rs  for  a  smgle  bottom 
bounce  only,  while  the  second  has  a  bottom  /  surface  /  bottom  bounce)  this  is  done  by 
pressing  the  Calculate  button.  The  positions  of  the  starts  and  ends  of  both  returns  ar 
found  by  searching  in  the  expected  regions  (see  appendix  5)  for  4  comecutrve  samples 
of  increasing  value  (a  backwards  search  is  performed  to  find  the  ends). 

Once  the  start  of  the  first  return  is  found  the  actual  seafloor  depth  is  calculated  and 
displayed  in  the  Seafloor  Depth  text  box.  This  new  seafloor  depth  is  then  used  m  the 
calculiion  of  the  depth  transformation.  Also  after  the  start  of  the  first  return  is  found, 
the  background  noise  level  is  removed.  This  is  done  by  averaging  over  the  ^e  befom 
the  first  return  and  subtracting  this  level  in  the  statistics  calralahons  The  pomts 
averaged  are  between  the  80*  sample  and  0.9*  the  start  time  of  the  first  return.  e 
average  is  not  taken  from  the  first  sample  due  to  a  ringing  from  the  echo  sounder 
appearing  in  the  first  few  samples. 


To  view  the  ping  returns  and  check  if  the  program  has  found  the  start  and  end 
positions  correctly  a  plot  of  the  ping  return  can  be  viewed  by  pressing  the  Graph  Input 
Wave  button  (see  Figure  6).  The  wave  plot  has  different  colours  showmg  which  rehirr^ 
were  found  in  the  return  signal.  The  red  line  is  the  original  signal,  the  green  Ime  is  the 
first  return  position  calculated,  and  the  blue  line  is  the  second  return  position 
calculated  (see  Figure  6).  The  green  and  blue  lines  should  encompass  the  red  hne  at  the 
return  sections.  If  no  return  is  found  or  a  section  of  it  is  missed,  then  the  part  ^ssed 
will  show  up  in  red  (see  Figure  7).  Also  there  is  a  vertical  blue  line  durmg  the  first 
return  to  show  the  start  of  its  tail  (see  Figure  6)  (the  tail  of  the  first  return  is  from  foe 
start  of  foe  first  return  +  foe  ping  duration  to  foe  end  of  foe  first  return  (see  Figure  9)). 
The  red  section  seen  at  foe  beginning  of  foe  plot  is  due  to  ringing  from  foe  echo 
sounder  ping,  it  is  not  a  bottom  return  (see  Figure  6).  The  horizontal  line  at  zero  (which 
is  normally  blue)  is  for  plotting  purposes  only  (see  Figure  6) . 
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Figure  6:  Returned  signal  for  the  actual  water  depth  with  both  returns  found  correctly.  This 
display  is  similar  to  what  is  seen  on  an  echo  sounder  display,  but  only  shews  one 
ping.  The  x  axis  is  depth  (m)  and  can  be  used  to  input  a  new  Seafloor  Depth  into 
the  Input  Values  frame  if  the  first  return  is  not  found  correctly.  The  y  axis  is 
power,  but  there  are  no  units  for  this  axis  since  the  signal  is  not  calibrated. 


Figure  7:  Returned  signal  when  the  start  of  the  first  return  is  not  found  correctly  and  the  second 
return  is  not  found  at  all.  The  red  line  at  around  12  m  shows  that  part  of  the  first 
return  was  missed.  The  red  section  at  30  m  shows  that  the  second  return  was 
missed. 
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If  a  return  is  not  found  because  of  noise  or  dropouts  then  the  signal  can  be 
smoothed  or  the  dropouts  removed  by  pressing  the  relevant  button.  Dropouts  should 
be  removed  before  the  signal  is  smoothed  since  smoothing  wiU  affect  the  dropout 
spikes.  The  smoothing  is  done  using  a  moving  3-point  filter  that  has  a  bias  of; 

Power(ii)  =  0.25  x  Power(ii-l)  +  0.5  x  Power(ii)  +  0.25  x  Power(ii+l)  (6) 

This  is  used  to  reduce  any  high  frequency  noise  in  the  sipal.  The  dropoute  are 
removed  by  looking  for  any  very  low  amplitude  signals  with  higher  amplitude  signals 

each  side. 

Also  if  the  seafloor  depth  has  changed  drastically  since  tiie  last  pmg  and  the  first 
return  can  not  be  found,  a  new  seafloor  depth  can  be  entered  in  the  Seafloor  Depth  text 
box.  This  enables  a  different  section  of  the  return  to  be  searched  for  the  start  of  the  nrst 
return,  (see  appendix  5) 

If  only  the  first  return  is  found,  statistics  for  the  first  return  are  calculated,  and  the 
FFT  and  Integrated  Power  plots  are  produced.  To  view  these  forms  just  press  on  the 
relevant  buttons. 

If  the  data  is  smoothed  too  many  times  the  original  waveform  can  be  recovered  by 
pressing  the  Refresh  The  Data  button. 

If  a  text  box  in  the  input  values  section  is  changed  the  current  data  and  any  future 
data  will  be  calculated  using  the  new  inputs.  The  Sounder  Depth  Setting  can't  be 
changed  since  the  echo  sounder  file  is  already  open.  If  this  was  entered  incorrectly  the 
program  has  to  be  restarted,  and  this  value  entered  before  the  echo  sounder  file  is 
reopened. 

After  this  ping  return  has  been  analysed  the  next  ping  can  be  loaded  by  pressing 
the  Sample  Next  Ping  button.  At  this  time  there  is  no  saving  of  the  statistical  data 
produced,  because  the  format  of  the  output  data  is  not  yet  determined. 
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3.3  Statistics  Produced 

First  Return  Statistics 

Once  the  start  and  end  points  have  been  found  for  the  first  return,  the  following 
calculations  are  done  on  the  signal  (see  Figure  10). 

For  the  whole  of  the  first  return; 

•  The  peak  power  and  peak  time. 

The  peak  power  and  times  are  found  by  looking  for  the  three  consecutive  points 
with  the  higjiest  sum,  between  the  start  and  end  points.  Times  are  calculated 
from  the  start  of  the  first  return. 

•  The  half  width. 

The  half  width  is  the  time  width  at  half  ttie  peak  power. 

•  The  centroid  position  as  power  and  time. 

The  centroid  Hmps  are  found  using  the  2nd  moments  (M.  R.  Spiegel  1972),  while 
the  centroid  powers  are  found  by  calculating  the  total  energy  of  the  return  and 
the  energy  from  the  peak  down  until  it  is  half  the  total  energy.  Times  are 
calculated  from  die  start  of  the  first  return. 

•  The  standard  deviation,  skewness,  and  kurtosis  of  power. 

The  standard  deviation.  Skewness,  and  Kurtosis  are  calculated  from  moments 
calculated  about  the  centroid  (M.  R.  Spiegel  1972). 
s2  =  m2 ,  a3=m3/  s^  ,  a4=m4/  s* 

where  s  =  Standard  deviation 

m2  =  2nd  moment  about  the  centroid 
m3  =  3rd  moment  about  the  centroid 
nu  =  4th  moment  about  the  centroid 
as  =  Skewness 
34  =  Kurtosis 

•  The  1st,  2nd,  3rd,  and  4th  moments  of  power,  calculated  from  the  start  of  the  first 

return. 

•  The  median  power  . 

•  An  Integrated  Power  Plot.  Power^  vs  Time  (sec)  (see  Figure  11). 

The  integrated  power  is  a  cumulative  sum  of  the  square  of  the  power 
normalised  to  one,  with  the  power  stunmed  from  the  start  of  the  return  time. 
(See  Lurton  and  Pouliquen  (1992)  for  examples.) 
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•  An  FFT  plot.  Power  VS  Frequency  (Hz)  (see  Figure  12). 

The  FFT  plot  was  calculated  using  a  Radix  2  decimation  (J.G.  Proakis  and  D.u. 

Manolakis  1989),  with  the  power  spectrum  (real2+imag2)  being  displayed. 

When  calculating  the  half  widths  and  centroid  powers,  the  width  of  the  return  at  a 
specific  power  level  has  to  be  found.  To  reduce  the  effect  of  the  noise  in  the  signal,  Ae 
curve  intercept  was  found  at  four  points,  two  points  each  side  of  the  return  s  peak  The 
inside  and  outside  values  on  each  side  were  foimd  and  then  averaged  (see  Figure  8). 


Figure  8 :  Intercept  point  calculation  for  one  side  of  the  peak. 


For  the  tail  of  tiie  first  return; 


•  The  centroid  power  and  time  (time  is  calculated  from  the  start  of  the  first  return). 


•  The  energy  of  the  taU.  CaUed  the  Tail  Ref.  1st  Return  for  the  reference  seafloor 

depth  and  Tail  1st  Return  for  the  actual  seafloor  depth,  in  the  energy  results 
frame.  The  energy  for  the  reference  seafloor  depth  can  be  compared  for 
different  sites.  The  actual  seafloor  depth  energy  can  not  be  used  for  comparison 
with  other  sites,  but  is  more  for  seeing  five  strength  of  the  signal  returned. 

The  tail  of  the  first  return  is  from  Ihe  start  of  the  first  return  +  the  ping  duration  to 
the  end  of  the  first  return  (see  Figure  9). 
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Figure  9:  Tail  Energy.  The  power  is  not  calibrated,  thus  has  no  units. 


3.3.1  Second  Return  Statistics 

For  the  whole  of  the  second  return; 

•  The  peak  power  and  peak  time  (time  is  calculated  from  the  start  of  tiie  second 

return). 

•  The  half  width  based  on  the  peak  power. 

•  The  centroid  position  as  power  and  time  (time  is  calculated  from  the  start  of  the 

second  return). 

•  The  energy  of  the  second  return.  Called  the  Ref.  2nd  Return  for  the  reference 

seafloor  depth  and  Actual  2nd  Return  for  the  actual  seafloor  depth,  in  the 
energy  results  frame. 

•  The  median  power. 
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Moments  for  the  First  Return  from  the  Start  of  the  Return — 
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Figure  10 :  Statistics  Form 
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3.4  Problems  foimd  from  Analysing  Test  Data 

The  digitised  test  data  is  mostly  at  zero  amplitude  when  an  echo  is  not  bemg  be^rd. 
This  shows  that  the  Echo  Listener™  is  not  recording  the  background  noise,  and  &us 
producing  an  artificially  higher  noise  floor.  The  cause  of  this  is  proMy  the  lack  of 
dynamic  range  when  the  depth  sounder  signal  is  being  digitised  ^e  current  eAo 
listener  has  a  12  bit  A/D  converter,  which  must  be  able  to  convert  the  rmgmg  at  the 
start  of  the  signal  without  an  over  load,  and  thus  is  not  recording  low  power  signals. 

The  start  and  end  points  of  the  second  return  are  currently  not  being  found  on  a 
high  percentage  of  pings.  This  is  due  to  the  low  power  levels  returned  to  the 
hydrophone  frL  the  second  return.  Since  the  Echo  Listener  «  currenfly  producmg  an 
artificial  high  noise  floor,  this  problem  could  cease  when  the  a^aal  noise  floor  is 
removed.  So  no  work  has  been  done  to  remove  this  problem  m  the  current  version  o 
ABCS,  and  this  should  be  looked  at  when  tiiere  is  a  better  signal  to  noise  ratio  for  the 

second  return. 


The  Echo  Listener  currently  being  used  was  designed  to  record  only  the  first  return, 
which  is  why  it  only  has  a  12  bit  A/D  converter.  This  can  be  used  to  do  acoustic  bottom 
classification,  since  the  QTC  System  uses  only  data  from  the  first  return  to  do  its 
classification. 


The  Echo  Listener  design  is  currently  beirig  updated  by  Sonar  Data  to  give  it  an 
extra  15  dB  dynamic  range.  Sonar  Data  have  tested  prototypes  of  the  new  desi^  and 
are  currently  having  the  new  circuit  boards  engineered  (Higginbottom  1999).  This  new 
model  should  remove  the  artificial  noise  floor  and  enable  the  second  return  to  be 

found. 
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4.  Conclusions 

The  oirrent  version  of  the  ABCS  program  has  manual  checking  for  the  start  and 
end  points  of  the  echo  sounder  returns.  The  first  returns  are  normaUy  found  correctly 
after  one  three  point  smoothing  of  the  raw  data,  but  are  sometimes  successfully  found 
with  no  smoothing.  The  second  return  is  normally  not  found,  due  to  a  low  signal  to 
noise  ratio  of  this  return,  caused  by  the  Echo  Listener™  currently  being  used  (see  3  4 
Problems  found  from  Analysing  Test  Data).  This  problem  can  either  be  fixed  or  the 
second  return  ignored.  The  first  return  has  enough  data  to  determme  the  acoustic 
bottom  classification,  since  QTCTm  uses  only  the  first  return  in  its  classification. 

When  the  digitised  echo  sounder  signal  has  a  higher  dynamic  range  or  it  is  decided 
to  ignore  the  second  return,  raw  data  should  be  obtained  at  a  number  of  sites  where  the 
sediment  type  is  known.  This  data  can  then  be  used  to  refine  the  ABCS  program,  so 
that  the  finding  of  the  returns  can  be  totally  automated,  without  the  need  to  check  if 
the  returns  were  foimd  correctly.  Also  a  format  for  the  statistical  data  produced  can  be 
decided  and  the  analysis  of  these  statistics  begun. 
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Appendix  1:  Returned  Ping  Shape  and  Energy 
Variations  with  Seafloor  Depth 

The  shape  and  energy  of  the  returned  signal  can  change  significantly  with  seafloor 
depth  (see  Figure  13)  as  a  result  of  ping  dilation  or  compression,  spreading  losses,  and 
absorption. 


Time  (sec) 


Figure  13:  Haw  seafloor  depth  effects  the  returned  signal  Calculated  using  Levels. 

Ping  Dilation  or  Compression 

When  the  depth  of  water  changes,  the  start  time  of  the  returned  ping  changes  along 
with  the  duration  of  the  ping.  The  changing  of  the  duration  of  the  returned  ping  is 
called  ping  dilation  or  compression,  depending  on  whether  the  ping's  duration  has 
increased  or  decreased. 
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To  see  how  tiie  returned  ping  can  change  duration  as  the  seafloor  depth  differs, 
consider  the  same  source  ping  insonifying  the  seafloor  at  two  different  depdis  to  an 
angle  of  20  (see  Figure  14).  The  wave  fronts  produced  by  the  source  will  be  spherical  in 
shape.  The  leading  edge  of  the  ping  will  be  insonifying  ffie  seafloor  at  an  angle  0and 
there  will  be  some  point  (point  A  in  Figure  14)  in  the  source  ping  that  is  just  about  to 
start  insonifying  the  seafloor,  with  time  between  these  two  points  being  St. 


Since 

dt-i  =  c  Sri  ,St2  =  c  Sri 


and 


Sr^ 


-Jl— 


d, 


X  ^2 


Cos9 

where  c  =  radial  speed  of  sound 

(S'  =the  distance  between  the  points 
d  =  depth 


eq(7) 


eq(8) 


Figure  14:  Two  different  depths  with  the  same  insonified  angle,  when  the  insonified  area  is  a 
circle. 


Combining  eq(7)  and  eq(8)  gives: 

d. 
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Therefore  the  time  taken  to  insonify  the  seafloor  to  a  specific  angle  after  the  start  of 
the  ping  first  hits  the  seafloor  takes  longer  for  deeper  water.  The  return  path  back  to 
the  hydrophone  will  have  a  similar  delay.  For  a  constant  sampling  rate,  more  samples 
will  be  taken  before  a  particular  angle  is  reached  for  a  deeper  bottom  compared  to  a 
shallower  bottom.  Therefore  the  deeper  the  water  die  longer  the  duration  of  the 
returned  pulse  (see  Figure  15).  This  result  also  holds  for  an  insonified  annulus  on  the 
seafloor. 


Figure  15:  Ping  variations  due  to  seafloor  depth  changes.  di  is  the  reference  seafloor  depth  and 
dz  is  the  new  seafloor  depth 

Spreading  Losses  and  Absorption 

An  echo  sounder  ping  experiences  only  significant  spherical  spreading  loss,  due  to 
the  short  distance  over  which  it  operates.  Since  spherical  spreading  is  dependent  on 
distance  the  amount  of  spreading  loss  will  vary  depending  on  depth. 

The  absorption  of  an  echo  sormder  ping  occurs  as  it  passes  through  the  water.  The 
amotmt  of  absorption  is  dependent  on  the  distance  travelled  through  the  water,  and 
the  temperature,  salinity,  and  pressure  of  die  water.  Therefore  as  the  seafloor  depth 
changes  die  amount  of  absorption  will  change,  but  since  the  amount  of  absorption 
occurring  in  the  distances  travelled  by  the  first  two  returns  of  an  echo  soimder  ping  is 
small3 ,  the  effect  of  absorption  is  minimal. 


3  For  a  temperature  of  20  ®C,  salinity  of  35,  and  depth  of  10  m,  die  absorption  at  10, 30,  50,  and 
100  kHz  is  0.761, 5.19, 13.0,  and  38.0  dB/km  respectively  (from  Francois  and  Garrison  1982) 
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Appendix  2:  Geometric  Errors  when  Converting  to  a 
Reference  Seafloor  Depth 


Consider  a  ping  of  duration  At  coining  from  an  omnidirectional  source.  This 
energy  will  radiate  out  from  the  source  hitting  the  seafloor,  and  some  of  the  energy  will 
then  be  backscattered  back  to  the  hydrophone.  Since  the  energy  emitted  from  the 
source  will  be  spherical  in  shape  and  only  At  in  duration,  some  of  the  beam  will  first 
insonify  a  circle  on  the  seafloor.  The  circle  will  widen  with  time  imtil  the  end  of  the 
ping  reaches  the  seafloor.  Then  an  insonified  annulus  will  form  on  the  seafloor,  which 
will  move  out  radially  from  the  centre. 

To  show  how  the  geometric  errors  occur,  consider  an  insonified  annulus  on  the 
seafloor  wifli  a  reference  seafloor  depth  of  50  m  and  an  actual  seafloor  depth  of  100  m 
which  is  transformed  to  the  reference  seafloor  depth. 

Using  a  realistic  case  where  the  ping  duration  At  =  2  msec,  $=  30°,  and  tiie  speed  of 
sound  c  =  1500  m/s  (see  Figure  21). 

For  the  100  m  case: 

0  /I 

Since  CosO  =  —  (eq(36)  appendix  4),  therefore  t=0.15396  sec  at  100  m. 
ct 

Also  Cos(d  -  A^)  = - -  (eq(37)  appendix  4) 

c{t  -  AO 

This  gives  A0=  1.33®  at  100  m 

d  t 

If  this  was  transformed  to  the  reference  seafloor  depth  of  50  m  using  t'=  — ^  the 

d 

time  would  become  t'  -  0.07698  sec,  but  0  and  A6  would  still  be  30°  and  1.33° 
respectfully.  This  gives  an  insonified  area  of  236  m*  at  the  reference  seafloor  depth. 

For  the  50  m  case: 

The  time  t  would  equal  0.07698  sec,  but  the  ping  duration  At  would  still  be  2  msec. 
Doing  the  above  calculation  again  gives  Ad=  2.76°  and  an  insonified  area  of  477  m^. 

This  shows  tiiat  the  size  of  the  insonified  annulus  on  the  seafloor  is  significantly 
different  when  comparing  the  100  m  seafloor  depth  adjusted  to  the  50  m  with  the 
reference  seafloor  at  50  m  depth. 
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Appendix  3:  The  'Levels'  Program 

The  Levels  program  was  written  to  calculate  the  expected  first  return  from  any 
shaped  source  ping  envelope.  The  return  could  then  be  transformed  to  a  reference 
depth  to  give  an  idea  of  how  well  a  specific  transform  is  working.  The  ping  envelope 
shape  used  to  simulate  the  Furuno  echo  sotmder  was  a  2  msec  half  sine  wave  envelope 
(see  Figure  2). 

Since  any  source  can  be  made  up  from  a  number  of  piece-wise  constant  power 
inputs  of  varying  amplitude  (see  Figtire  2),  and  the  return  from  a  constant  power 
source  can  be  calculated  using  just  two  equations  (see  appendix  4  for  the  derivation  of 
die  two  equations  below),  the  program  was  written  to  calculate  the  return  using  this 
method. 


(eq(28)  and  eq(31),  appendix  4) 


(eq(38)  and  eq(39),  appendix  4) 


t  =  Time  since  the  start  of  the  constant  power  input  was  produced  by  the 
source. 

P  =  Power  coming  from  the  omnidirectional  source. 

R  =  Power  received  by  the  hydrophone. 

=  The  backscatter  ratio. 

At  -  The  duration  of  the  constant  power  input. 

The  levels  program  uses  the  above  two  equations  to  calculate  the  return  for  one 
constant  power  input  of  duration  At  (see  Figure  16).  It  then  simulates  the  return  for  the 
whole  ping  by  adjusting  the  amplitude  and  doing  a  time  shift  of  At  for  each  of  the 
constant  power  inputs  which  make  up  the  ping  (see  Figtire  2).  All  the  inputs  are  then 
added  to  give  the  returned  signal  (to  give  a  better  idea  of  how  the  inputs  are  added  the 
first  four  input  returns  are  shown  in  Figure  17). 


R  = 


The  two  equations  used  were: 

PP 


6d^ 


1- 


ct  J 


V  V  -  -  y 

,  ,  .  .  ,  2d  .  ^2d  . 

which  IS  used  when  — <t<  —  +  m 
c  c 

And 


R  = 


APpd 


f 


3c" 


1 


1 


{t-Aif  f 


which  is  used  when  t  ^ 


2d 


■At 


where  d  =  Depth  to  the  seafloor. 

c  =  Sound  speed  in  the  water. 
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Figure  17:  Summing  the  first  four  constant  power  inputs  for  the  sine  wave  used  (see  Figure  2). 

There  are  more  inputs  to  simulate  the  whole  sine  input,  hut  these  were  left  out  so 
the  graph  is  not  too  complicated.  Each  input  had  a  duration  At  =  0.01  msec.  Input 

1  had  zero  amplitude.  Input  2  had  a  time  shift  of  At.  Input  3  had  a  time  shift  of 

2  At.  Input  4  had  a  time  shift  of  3 At. 

When  calcnilating  the  returns  the  following  assumptions  were  made: 

•  The  seafloor  is  flat. 

•  The  backscatter  strength  is  constant  with  angle. 

•  An  omnidirectional  source  and  receiver. 

Thus  any  plots  produced  by  this  program  will  have  these  assumptions. 
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Appendix  4:  Calculation  of  the  Power  Returned  from  a 

Constant  Power  Input 


Consider  a  constant  power  ping  of  duration  At  coming  from  an  omnidirectional 
source.  This  energy  will  radiate  out  from  the  source  hitting  the  seafloor,  and  some  of 
the  energy  will  then  be  backscattered  back  to  the  hydrophone.  Since  the  energy  emitted 
from  the  source  wiU  be  spherical  in  shape  and  only  At  in  duration,  some  of  the  beam 
will  first  insonify  a  circle  on  the  seafloor.  The  circle  wiU  widen  with  time  until  the  end 
of  the  ping  reaches  the  seafloor.  Then  an  insonified  annulus  will  form  on  the  seafloor, 
which  will  move  out  radially  from  the  centre. 

Calculation  of  the  Power  in  a  Beam  from  an  Omnidirectional  Source 

Consider  the  energy  emitted  from  the  omnidirectional  source  of  constant  power 
and  a  ping  duration  of  At.  This  energy  will  radiate  away  from  the  source  forming  a 
spherical  shell  of  width  At  and  radius  b  (assuming  a  constant  sound  speed  in  the 
water).  The  power  level  emitted  by  the  source  will  be  equal  to  the  total  power  over  the 
surface  area  of  the  sphere.  Since  the  source  is  omnidirectional  the  power  wiU  be  evenly 
distributed  over  the  sphere. 

eq(10) 

Power  level  per  unit  area  over  the  sphere. 

Radius  of  the  sphere. 

Power  level  emitted  from  the  source. 

Therefore  the  power  level  m  a  beam  of  half  width  ^  (see  Figure  18)  will  be: 

B  =  P^a  eq(ll) 

where  B  =  Power  level  m  the  beam. 

cr=  Surface  area  of  the  spherical  beam  at  radius  b,  with  half  width 


where  Ps- 
&  = 
P  = 
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Figure  18:  Beam  of  half  width  ^and  duration  At. 

The  equation  describing  the  spherical  beam  is: 

-x^  -y^  eq(12)  (see  Figure  18) 

The  surface  area  (cr)  of  the  spherical  beam  of  half  width  ^  will  be: 

C7=  jj-Jl  +  f^+fydA  eq(13)  (Grossman  1984) 

where  =  The  region  defining  the  beam  at  radius  b 

j  y  r  2  2  2 

^  b  -X  -y 

Therefore 

i^P^-x^-y^ 
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Substituting  y=nSina>  and  x=fiCosco  gives: 


nllt 


rdcodji 


eq(15) 


where  ju  =  Distance  from  frie  z  axis 

(O  =  Angle  from  the  x  axis,  centred  on  the  z  axis 
ie.  a  cylindrical  coordinate  system 


Solving  this  gives: 


eq(16) 


Since  a  =  bSin^  (see  Figure  18)  therefore: 

O  =  2jtb^i\-Cos^)  eq(17) 


Combining  eq(10),  eq(ll),  and  eq(17)  gives: 

B  =  ^{i-Cos^)  eq(18) 

where  P  =  The  power  level  emitted  by  the  source. 

B  =  The  power  level  of  a  beam  of  half  width 

Calculating  the  Insonified  Circle  Power  Backscattered. 

An  insonified  circle  on  the  seafloor  will  exist  from  the  time  tfie  front  of  the  ping 
first  reaches  the  seafloor,  imtil  the  end  of  the  ping  reaches  the  seafloor. 

d  d 

ie.  when  —<t<  — h  At  eq(19) 

c  c 

where  d  =  Depth  to  the  seafloor. 

c  =  Soimd  speed  in  the  water. 

t  =  Time  since  the  start  of  the  ping  was  produced  by  the  source. 

At  =  The  duration  of  the  ping, 
assuming  a  constant  sound  speed  and  a  flat  seafloor. 
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Figure  19:  Insonified  circle  on  the  Seafloor 

The  level  of  insonification  in  the  circle  (shown  in  Figure  19)  will  be  equal  to  the 
power  in  the  angular  beam  of  width  2^  which  produced  the  circle. 

Therefore  the  level  of  insonification  in  the  circle  will  be: 

B  =  Cos  0)P  eq(20)  (from eq(18)) 

where  P  =  Power  level  of  the  ping  coming  from  the  source. 

B  =  Level  of  insonification  in  the  circle,  (see  Figure  19) 

Since  a  spherical  coordinate  system  was  used  to  calculate  eq(20),  spherical 

The  size  of  the  insonified  circle  increases  with  time  (ie.  ^  will  increase  with  time), 
but  only  exists  between  the  two  constraining  times  (eq  (19)).  The  level  of  insonification 
will  vary  over  the  circle,  with  the  strongest  insonification  in  the  centre  of  the  circle. 

Since  eq(20)  is  the  level  of  insonification  in  the  circle,  and  the  circle  is  increasing  in 

diameter  with  time,  —  will  equal  the  level  of  insonification  at  the  outside  edge  of  the 

dO 

insonified  circle  (this  is  only  true  for  a  constant  power  ping,  since  it  requires  that  any 
point  in  the  insoiiified  circle  have  a  constant  insonification  level  with  time). 

*  For  a  constant  angle  beam,  the  amount  of  energy  in  the  beam  will  be  constant  if  there  is  no 
absorption  in  the  water.  This  is  because  the  diameter  of  the  beam  increases  with  range  at  the 
same  rate  as  the  signal  at  any  one  point  decreases. 
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R^^Ei.\sm<l,.d(l>  eq(23) 

where  R  =  Power  level  received  at  tiie  hydrophone. 


Since  eq(24)  (see  Figure  20) 

d 

Therefore 

R=  ^^^Sin^.Cos^^.d<l)  eq(25) 

Assuming  the  backscatter  ratio  is  constant  with  incident  angle  gives; 


R  =  ^  Sin^.Cos^^.d<l> 


Therefore 


eq(26) 


R^lLi^-Cos^e)  eq(27) 


This  signal  would  arrive  at  the  hydrophone  between  the  times; 

eq(28) 

c  c 

d 

Since  Cos9  =  — 
r 

eq(29) 

and  r  =  — 

2 

eq(30) 

where  r  =  The  distance  travelled  by  the  start  of  the  ping  to  the  seafloor  (see  Figure  19). 
t  =  The  travel  time  to  the  seafloor  and  back  to  the  hydrophone. 


Therefore  the  power  level  received  by  the  hydrophone  at  any  instant  during  the 
times  specified  by  eq(28)  will  be; 


R 


P£ 

6d^ 


1- 


ct 


where  t  =  Time  the  backscattered  signal  is 


eq(31) 

received  by  the  hydrophone. 
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The  power  backscattered  by  the  insonified  annulus  (see  Figure  21)  to  the 
hydrophone  at  a  time  given  by  eq(32)  will  be; 

Sin^.Cos'iff.d^  eq(33)  (see  eq(26)) 

2d^ 

R  =  -^((l  -  Cos  ^0)-  (1  -  Cos  ^{0  -  A0)))  eq(34) 

6d 

Therefore 


R  = 


E£^{CosH0-A0)-Cos^0) 


eq(35) 


Since 

ct 

”T 


Cos0  =  — 


c{t  -  At) 


and  Cos{0  -  A0)  =  — 


Therefore 

Cos0  =  — 
ct 

and  Cos{0  —  A0)  =  — - —  6q(37) 

c{t  -  At) 

Combining  eq(35),  eq(36),  and  eq(37)  gives  the  power  level  returned  to  the 
hydrophone  for  an  annulus. 


1 _ ll 

3c"  [{t-Atf  t^J 
2d 

which  occurs  when  t  ^ - H  At 

c 


eq(38) 

eq(39) 
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Appendix  5:  Search  Areas  used  by  'ABCS'  to  find  the 

Seafloor  Returns 

The  ABCS  program  uses  the  seafloor  depth  to  calculate  the  expected  regioiis  in 
which  the  start  and  finish  of  the  returns  should  be  (see  Figure  22).  To  find  a  start  or  end 
point  the  program  looks  for  4  consecutive  points  of  increasing  value  in  the  expected 
regions. 


The  first  point  searched  for  is  the  start  of  the  first  return.  This  search  starts  from 
0.7*Start  time  estimated  from  the  seafloor  deptti  and  goes  along  a  time  interval  of  the 
sample  duration/ 5.0.  The  time  estimated  from  the  seafloor  depth  =  2*seafloor 
dep&/soimd  speed.  If  the  seafloor  deptii  in  this  text  box  is  off  by  a  large  amoxmt  the 
expected  region  might  not  contain  the  start  point.  If  a  start  point  is  found,  the  seafloor 
depth  is  calculated  from  this,  with  all  the  other  search  regions  using  this  new  seafloor 
depth.  If  this  start  point  is  not  fotmd  no  other  points  are  looked  for. 

The  next  point  searched  for  is  the  start  of  the  second  return.  This  search  starts  at 
1.9*Start  of  the  first  return  and  has  a  seeirch  duration  of  the  sample  duration  /  5.0. 
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The  end  of  the  returns  are  looked  for  next,  and  are  searched  for  in  the  reverse 
order.  With  these  points  the  times  searched  depends  on  whether  the  second  start  point 
was  found.  If  the  second  start  point  was  found,  the  end  of  the  first  return  is  searched 
for  between  the  two  start  points,  and  the  end  of  the  second  return  is  searched  from  the 
end  of  the  samples  for  the  ping  to  die  start  of  the  second  return.  If  the  start  of  the 
second  return  is  not  foimd,  the  end  of  the  first  return  is  searched  for  from  the  end  of 
die  ping  samples  to  the  start  of  the  first  return,  with  no  search  for  the  second  return 
end. 

The  background  noise  level  is  also  calculated  using  the  start  time,  this  value  is 
taken  as  an  average  of  the  points  between  the  80*  sample  and  0.9*  the  start  time  of  the 
first  return.  The  average  is  not  taken  from  the  first  sample  due  to  a  ringing  in  the  echo 
sounder  appearing  in  the  first  few  samples. 
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